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1 The problem 

Air Traffic Flow Management is the regulation of air traffic in order to avoid exceed- 
ing airport or flight sector capacity in handling traffic, and to ensure that available 
capacity is used efficiently. The number of flights taking off or landing from a certain 
airport or the number of planes traveling in a particular sector are functions of several 
variables including: 

1. The number of runways available 

2. ATC capacity 

The Air Traffic Controllers and the Control centre perform various tasks like: 
Terminal Approach, Final Approach, Ground (arrivals), Gate operations (ar- 
rivals), Airport Surface, Gate operations (departures), Ground (departures) and 
finally the Take-Off / Transition to Terminal and Center. Thus the number and 
expertise of ATCs available determines the airport capacity to a certain extent. 

3. Restrictions as to which aircraft can follow an aircraft of a given class 

This is because an aircraft periodically sheds vortices from the wing tips. These 
vortices can interact with the boundary layer around the wing of an aircraft 
following in the first aircraft's wake and destroy the lift. The effect is more 
pronounced when the following aircraft is of a smaller size class then the leading 
one. This means that a delay must be imposed such that F t > L t + 5f,l where 
Ft and L t denote the times at which the follower and the leader take-off or land 
and Sf,l denotes the delay. 

4. Airspace restrictions 

Aircraft have to follow certain "corridors" when traveling in air. This restricts 
the number of aircraft that can be airborne at any given instant. Finally, most 
of the restrictions mentioned above, are functions of time too. This makes the 
problem dynamic in nature. 

2 Existing models 

Airport and airspace capacity is the major cause of congestion. The difficulty in 
dealing with this parameter is its uncertainty heavily influenced by weather conditions 
among other factors. Short-term solutions (12 hours time horizon) for air traffic flow 
management include ground-holding policies. These policies are motivated by the 
fundamental fact that airborne delays are much costlier than ground delays, since the 
former include fuel, maintenance, depreciation and safety costs. Therefore, the aim 
of ground-holding policies is to translate anticipated airborne delays to the ground. 
In addition to ground holding, distributing the traffic efficiently across the airspace 
helps in reducing congestion in the network. 
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Several models have been proposed in the literature for problem solving under 
different real-life hypotheses. The Single- Airport Ground Problem (SAGHP) only 
considers a single airport and the goal is to produce ground-holding schedules. The 
Multi- Airport Ground-Holding Problem (MAGHP), Andreatta et.al.[l] considers a 
network of airports such that the ground-holding policies for one of them have im- 
pact on the other airport schedules, as can be seen in Vranas et.al. [ I] . The Air 
Traffic Flow Management Problem (TFMP) is similar to MAGHP but also considers 
the airspace network besides the airport capacity, Bertsimas et.al. [2]. The main 
assumption on the problem made by most of the approaches described in the open 
literature is the deterministic character of the data. 

3 Motivation and Objectives 

The present work focuses on the seminal paper, due to Bertsimas et.al [2] in which 
the polyhedral structure of the TFMP is analyzed. The authors of this paper claim 
that by using a transformation of the decision variables, a formulation in which some 
inequalities define the facets of the convex hull of the TFMP was obtained. 

This is a significant improvement in the formulation, since Polyhedral theory 
indicates that under the conditions mentioned above, the LP relaxation of the TFMP, 
might be integral. This means that we may not need techniques like branch and 
bound, at least in some of the cases. In this work an attempt has been made to 
understand and verify the claims made in the paper [ ]. To be specific: 

1. The physical problem was studied in terms of the several formulations 
(SAGHP, MAGHP, TFMP) available. 

2. We tried to understand the main results of of Polyhedral Combinatorics, prin- 
cipally through the textbook by Wolsey and Nemhauser [3]. 

3. The paper due to Bertsimas et.al [2] was examined in the light of this theory 
to understand how the formulation helps to achieve integral solutions. 

4 Method 

4.1 The Traffic Flow Management Problem (TFMP) 

This subsection has been reproduced verbatim from Bertsimas et.al [2]. Consider a 
set of flights, & = {1, a set of airports, = {1, ...,K}, a set of time peri- 

ods, = {1,...,T}, and a set of pairs of flights that are continued, ^ = {(/,/) : 
/ is continued by flight/}. We shall refer to any particular time period t as the "time 
f ' . The problem input data are given as follows: 
Nf = number of sectors in flight f's path 
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P(f, i) = the i th sector in flight f's path 

P f = {P(f,i):l<i<N f ] 

D k (t)= departure capacity of airport k at time t 

Ak(t)= arrival capacity of airport k at time t 

Sj(t) = capacity of sector j at time t 

df = scheduled departure time of flight / 

rf — scheduled arrival time of flight / 

Sf = turnaround time of an airplane after flight / 

= cost of holding flight / on the ground for one unit of time 
c[ = cost of holding flight / in the air for one unit of time 
/ fj = number of time units that flight / must spend in sector j 
Tj = set of feasible times for flight / to be in sector j 

Note that by "flight," we mean a "flight leg" between two airports. Also, flights 
referred to as "continued" are those flights whose aircraft is scheduled to perform a 
later flight within some time interval of its scheduled arrival. 

Objective: The objective in the TFMP is to decide how much each flight is going 
to be held on the ground and in the air in order to minimize the total delay cost. We 
model the problem as follows. 
Decision variables: 

j J 1 if flight / arrives at sector j by time t 
ft \ otherwise 

Note that the wj t are defined as being 1 if flight / arrives at sector j by time t. 
This definition using by and not at is critical to the understanding of the formulation. 
Also recall that we have also defined for each flight a list Pf of sectors which includes 
the departure and arrival airports, so that the variable w^ t will only be defined for 
those sectors j in the list Pf. Moreover, we have defined Tj as the set of feasible 
times for flight / to be in sector j, so that the variable Wj t will only be defined for 
those times within Tj . Thus, in the formulation whenever the variable wi t is used, it 
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is assumed that this is a feasible (f,j,t) combination. Furthermore, one variable per 
flight-sector pair can be eliminated from the formulation by setting w 3 - j = 1 where 

fT f 

Tj is the last time period in the set Tj . Since flight / has to arrive at sector j by the 
last possible time in its time window, we can simply set it equal to one as a parameter 
before solving the problem. 

Having defined the variables we can express several quantities of interest as linear 
functions of these variables as follows: 

1. Noticing that the first sector for every flight represents the departing airport, 
then the total number of time units that flight / is held on the ground is the 
actual departure time minus the scheduled departure time, i.e., 

9f= Yl - w h-i) - d f 

teT*,k=P(f,i) 



2. Noticing that the last sector for every flight represents the destination airport, 
the total number of time units that flight / is held in the air can be expressed 
as the actual arrival time minus the scheduled arrival time minus the amount 
of time that the flight has been held on the ground, i.e., 

a f= Yl ^ w ft - w U-i) ~ r f-9f 

t€T*,k=P(f,N f ) 

The objective of the formulation is to minimize total delay cost, and the TFMP is 
hence: 

TFMP : Min z = <%9f + c } a f 
subject to Y ( w ft - w ),t-i) < Dk{t),Vk eJtT,te ST (1) 

/:P(/,l)=fc) 

J2 {w k ft -wl t ^)<A k {t)^ke,je,te3r ( 2 ) 

f:P(f,N f )=k) 

(vb-vlJzSMMeSitef (3) 

f:P(f,i)=j,P(f,i+l)=j 

j < <n lVfe*,teT},j = P(f,i), 



^^^°'ii = m ? + i),<^ (4) 



WU ""'I P(f,l) = k = P(f,N f 



(5) 
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<.,■-<., , > 0,V/ e eP f ,te& (6) 
w j ft e{o,i},Vfe&,jeP f ,te& (7) 

The first three constraints take into account the capacities of various aspects of 
the system. The constraint (1) ensures that the number of flights which may take 
off from airport k at time t, will not exceed the departure capacity of airport k at 
time t. Likewise, the constraint (2) ensures that the number of flights which may 
arrive at airport k at time t, will not exceed the arrival capacity of airport k at time 
t. In each case, the difference will be equal to one only when the first term is one 
and the second term is zero. Thus, the differences capture the time at which a flight 
uses a given airport. The constraint (3) ensures that the sum of all flights which may 
feasibly be in sector j at time t will not exceed the capacity of sector j at time t. This 
difference gives the flights which are in sector j at time t, since the first term will be 
1 if flight / has arrived in sector j by time t and the second term will be 1 if flight 
/ has arrived at the next sector by time t. So the only flights which will contribute 
a value 1 to this sum are the flights that have arrived at j and not yet departed by 
time t. Constraints (4) represent connectivity between sectors. They stipulate that 
if a flight arrives at sector j by time t + Ifj, then it must have arrived at sector j by 
time t where j and j are contiguous sectors in flight f's path. In other words, a flight 
cannot enter the next sector on its path until it has spent i^-time units (the minimum 
possible) traveling through sector j, the current sector in its path. Constraints (5) 
represent connectivity between airports. They handle the cases in which a flight is 
continued, i.e., the flight's aircraft is scheduled to perform a later flight within some 
time interval. We will call the first flight / and the following flight /. Constraints 
(5) state that if flight / departs from airport k by time t, then flight / must have 
arrived at airport k by time t — sj. The turnaround time, sy, takes into account the 
time that is needed to clean, refuel, unload and load and further prepare the aircraft 
for the next flight. In other words, flight / cannot depart from airport k, until flight 
/ has arrived and spent at least Sj time units at airport k. Constraints (6) represent 
connectivity in time. Thus, if a flight has arrived by time t, then Wft, has to have a 
value of 1 for all later time periods, t >t. 

4.1.1 Complexity of the TFMP 

Bertsimas et.al., [2], proved that the TFMP with all capacities equal to 1 is NP-hard, 
by showing equivalence with the Job-shop scheduling problem. 

4.1.2 Computational results due to Bertsimas et.al. [2] 

Bertsimas et.al, solved many instances of the MAGHP and the TFMP using a CPLEX 
2.1 solver. They used as many as 1000 flights, varying the number of connected flights 
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from 1/5 to 4/5 of the total flights being flown. The problem was solved over a 24-hour 
period, with intervals of 15 min. Even at the infeasibility border (which represents 
the data set such that even a minor change in parameters will make the constraints 
inconsistent), only 4% of the solutions were found to be non-integral, that too for the 
MAGHP. All instances of the TFMP had integral solutions. 

5 Polyhedral combinatorics to identify Facet defin- 
ing constraints 

The use of the variables denoting whether a flight has arrived by a certain time rather 
than at a certain time has advantages, achieved in terms of a "stronger" formulation. 
Key results from Polyhedral theory, relevant to the present work are included here and 
are used in sketching a constructive proof of the high incidence of integral solutions 
even for the LP relaxation of TFMP. 

5.1 Basic Concepts from Polyhedral Theory 

The definitions and theorems in this section, closely follow the development given in 
Wolsey and Nemhauser [3]. 

Definition 1. A set of k vectors, G $l n is defined to be affinely independent if 

k k 

\iX { = 0, = A, = Vz = 1, k 
i=i i=i 

Definition 2 (Polyhedron). A polyhedron P C dt n is a set of points that satisfy a 
finite number of linear inequalities; that is, P = {x G 3? n : Ax < b}, where (A, b) is 
an m x (n + 1) matrix. 

Definition 3. The dimension of a Polyhedron P is one less than the cardinality of 
the maximal set of affinely independent points in P. 

Lemma 1 (Dimension of P). If P C 3? n , then dim(P) + rank(A = ,b = ) = n where 
(A = ,b = ) is the equality set consisting of the rows corresponding to a l x — bi Vi G P. 

Definition 4 (Valid Inequality). An inequality ttx < ttq is called valid for polyhe- 
dron P C 3? ra if it is satisfied by all points of P. 

Definition 5 (Face). If ttx < tc q is a valid inequality for P and F = {x G P : nx = 
7To} then F is called a face of P, represented by (tt, 7To). (F ^ iff max{7c x : x G 

P} = TTo). 

Definition 6 (Facet). A face of P, is called a facet if dim(F) = dim(P) — 1. 
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Property 1 (Description of P). The facets of P are necessary and sufficient for 
the description of P. 

Definition 7 (Convex Hull). The set S is defined as S = P f] Z n where P C 3? n . 
The convex hull of S is denoted as conv(S) and conv(S) = {y G 3? n : y = £) i=1 ^i^*}; 
where i' 6 5 and A, G 3? + Vi, with 52 i=1 Aj = 1. 

Convex hull of IP 

— " — 1 1 - Feasible region, P 













s 










s, 

\ 








V \ 


V 



12 3 4 



Figure 1: The dash-dotted line indicates the feasible set P, while the solid line denotes 
the convex hull of the corresponding integral polyhedron on 3? 2 . An extremal point 
y G conv(S) has the property that $ x l ,x 2 G conv(S) such that y = \x x + \x 2 . 

Property 2. If n x < 7r defines a face of dimension k — 1 of conv(S) , 3 k affinely 
independent points x 1 , x k G S such that nx l = n , for i = 1, k. 

Definition 8 (Extreme point). A point x G P is called an extreme point if $x 1 , x 2 G 
P such that x = \x x + \x 2 . 

Property 3. A maximal valid inequality for S, is one that dominates all other valid 
inequalities. The set of maximal valid inequalities of S contains all the facet- defining 
inequalities of conv(S) , p. 207 [ }. 

The development in Wolsey [3] , uses a theorem relating the extremal points of the 
convex hull of S, to those of the set S itself. They state a property without proof, for 
which the following intuitive proof can be given. The point y G conv(S) is assumed 
to have the property that $x 1 ,x 2 G conv(S) such that y = ^x 1 + \x 2 . There are 
several cases depending on where the point y lies in Figure 1: 

1. y G conv(S')\{boundary of conv(S)}. This contradicts the assumption that y 
is an extremal point of conv(S). This is because all strictly interior points can 
be obtained as (non-trivial) convex combinations of some other points. 
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2. y G {boundary of conv(S)}\{x : x is an extremal point ofS"}. Refer to Fig- 
ure 1. However, such points on the boundary of conv(S) can be obtained 
as (non-trivial) convex combinations of the extremal points of 5*, (and S C 
conv(S)), thus contradicting the assumption. 

If y does not lie in any of the regions mentioned above, then y G S and S C 
conv(S). Now, if y is not an extremal point of S, 3 x 1 , x 2 , such that y = \x x + \x 2 , 
where x l ,x 2 G conv(S), contradicting the assumption again leading to the following 
lemma. 

Lemma 2. Hence, if such an extremal point of conv(S) 3, it must be one of the 
extremal points of S . 

Theorem 1. If P = {x G 9ft n : Ax < b} is the feasible set of the LP relaxation 
of the IP max{cx : x G S} where S = P fl Z n and the LP on the convex hull of 
S, max{cx : x G conv(S)} has an optimal solution, x* (which will be an extremal 
solution, since conv(S) is a polyhedron) , then x* will be an optimal solution to the IP 
on S too. 

5.2 Why the formulation used by Bertsimas et.al., [ ] is strong? 

We shall now try to explain the logic behind the occurrence of integral solutions 
when using the particularly strong formulation, given by Bertsimas et.al, [2] using a 
constructive proof. 

5.2.1 Constructive Proof 

(i) Try to find the conv(S), which is difficult usually. If we can find the conv(S) 
uniquely, we are done by Theorem 1. Else . . . 

(ii) Check which of the inequalities governing the IP are facets of the IP. This can 
be done by choosing each a % x < 6j in turn and finding the dim(F) where F is 
the set of points in IP satisfying a l x < bi as equality. If dim(F) = dim(IP) — 1, 
F is a facet according to Definition 6. 

(iii) The facets of S are also facets of conv(S) by Property 3. If, by using this 
iterative method, we can find the complete set of facets of conv(S), we are 
done, as solving the LP on conv(S) for an optimal solution (if it exists) will 
give an optimal solution to IP too, by Theorem 1. However, even if we cannot 
get all the facets, of conv(S) as is the case in real-life instances, the formulation 
is strong because of some constraints being facets and there is a high possibility 
of getting integral solutions. 
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6 Implementation 

We implemented the TFMP on artificially constructed flight data sets, to assess the 
computational performance of the model. Gnu Mathprog was used as the translator 
as it employs symbolic algebraic notation, which has advantages over languages like 
LINDO, which require that each constraint, alongwith the numerical values of the 
parameters be explicitly written down. The model was solved using the Open-source 
GLPK [6] linear solver. The standalone solver contained within GLPK, glpsol was 
used for this project. 

6.1 The Gnu Mathprog code 

Gnu Mathprog is an open-source code based on the C language. It is in fact the old 
version of AMPL as given in the paper by Kernighan et.al. [5]. The data consists of 
three primary sets: 

1. The set of airports and sectors, denoted as ' . 

K = {Mumbai, Bangalore, Pune, Goa, Calicut, Belgaum,Coimbatore} 



2. The set of flights, denoted as & . 

K = {Mu_Pu_Go, Go_Ca, Pu_Be_Ba, Go_Co_Ba } 

The Mu_Pu_Go flight is continued by the Go_Ca flight. Other flights are inde- 
pendent of each other. 

3. The set of time intervals, ST over which the schedule is to be prepared. 

The Gnu Mathprog model for the TFMP is shown below. The model is mostly 
self-explanatory, but the reader is referred to Kernighan et.al [5] for details. 

6.2 Gnu Mathprog Model 

/* 

GMPL model for the TFMP. 
Airports are: 

Mumbai, Bangalore, Pune, Goa, Calicut, Belgaum, Coimbat 
ore 

Flights are: 

* Mumbai-Pune-Goa connected to Goa-Calicut 
*Pune-Belgaum-Bangalore 
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Figure 2: The set of 7 flight sectors and 4 flights for South Western India. The 
Mumbai-Pune-Goa flight is continued by the Goa-Calicut flight. There are two other 
independent flights, Goa-Coimbatore-Bangalore and Pune-Belgaum-Bangalore. This 
is a hypothetical set of flights designed to illustrate concepts. 
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*Goa-Coimbatore-Bangalore 
Time slots are t=0 to 10. 

The chosen data set illustrates how the Pune-Belgaum-B 
angalore fit. 

is delayed because of restricted landing capacity at Ban 
galore airport 



*/ 



set N; 

#set of natural numbers 
set K; 

/*Airports or sectors*/ 
set F; 
/*Flights*/ 
set T; 

/*Time periods*/ 

param arr{k in K, t in T}; 
/*arrival capacity*/ 

param cap{k in K, t in T}; 
/* capacity*/ 

param atleast_one{f in F, k in K, t in T}; 
#flt. MUST reach sector k by time t 

param dep{k in K, t in T}; 
/*departure capacity*/ 

param sched_dep{f in F}; 
/*scheduled dep*/ 

param sched_arr{f in F}; 
/*sched arr time*/ 
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/*Sector flight transit minimal time*/ 
param tr{f in F,k in K} ; 

param turn{f in F}; 

/*turnaround time for an aircraft*/ 

param ground-[f in F}; 

/*ground holding cost per unit of time*/ 
param air{f in F}; 

/*air holding cost per unit of time*/ 

param feas{t in T, k in K, f in F}; 

#feasible flight-sector-time combinations 
param connect{g in F, f in F}; 
#set of connecting flights 
param sec_ord {f in F , n in N} symbolic ; 
#represents the n th sector in flight f 's path 

param maxi{f in F}; 

#cardinality of the set of sectors for a flight 
#####VARIABLES###### 

var w{f in F, t in T , k in K}, >=0, <=1; 
#=1 if fit. f arrives at sector k BY time t 
var grdel_fl{f in F} ; 
#ground delay for flight f 
var airdel_fl{f in F} ; 
#air delay for flight f 

#air delay for flight f 

var cost;# cost variable =max(0, actual cost) 
var dummyob j ; #dummy objective function 

minimize dummy : dummyob j ; #the objective function 

s.t. dummyobjl: dummyob j =sum{f in F} (ground [f] *grdel_fl [f] +air [f] *airdel_fl [f] ) ; 

s.t. pos:cost>=0; 
s.t. posl : cost>=dummyobj ; 
# the objective function 
s.t. grdel{f in F}: 
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grdel.f 1 [f ] =sum{t in T : f eas [t , sec.ord [f , 1] , f ] =1 

}t*(w[f ,t ,sec_ord[f , 1] ] -w[f ,t-l,sec_ord[f , 1] ] ) -sched_dep [f ] ; 

#ground delay for flight f 

s.t. ardel{f in F}: 

airdel_f 1 [f ] =sum{t in T : f eas [t , sec_ord [f , maxi [f]] ,f]=l } 
t* (w[f , t , sec_ord [f ,maxi [f ] ] ] -w [f , t-1 , sec_ord [f , maxi [f ] ] ] ) 
-sched_arr [f ] -grdel_f 1 [f ] ; 
#air delay for flight f 

s.t. domain{f in F,t in T, k in K : f eas [t ,k,f ] =0 . 0} : 
w[f ,t,k] =0; 

# no arrival/dep at infeasible combinations 

s.t. atleast{f in F, k in K, t in T: atleast_one [f ,k,t] =1} : 
w[f ,t,k]=l; 

# cap on time at which flights must arrive at a sector 

subject to depart {t in T, k in K } : 
sum{f in F : k=sec_ord [f , 1] and f eas [t ,k,f ] =1 } 
(w[f ,t,k]-w[f ,t-l,k] )<= dep[k,t] ; 

#departure capacity of airport 
s.t. arrive {t in T,k in K} : 
sum{f in F : f eas [t ,k,f ] =1 and k=sec_ord[f ,maxi [f ] ] } 
(w[f ,t,k] -w[f ,t-l,k] )<=arr [k,t] ; 
#arrival capacity of airport 
/*Sector capacity*/ 

s.t. capa {t in T,k in K} : sum{f in F, n in N :n<maxi[f] and 
f eas [t ,k,f ] =1 and sec_ord [f ,n] =k } 

(w [f ,t , sec_ord [f ,n] ] -w [f ,t , sec_ord[f ,n+l] ] )<=cap [k, t] ; 
/*Sector connectivity*/ 

s.t. seccon {f in F,n in N, t in T: n<maxi[f] 
and f eas [t , sec_ord [f ,n] ,f ] =1 } : 
(w [f , t+tr [f , sec_ord [f , n] ] , sec_ord [f , n+1] ] 
-w [f ,t , sec_ord [f ,n] ] )<=0; 

s.t. turnaround {f in F, g in F,t in T: 
f eas [t , sec_ord [f , 1] , f ] =1 and 
connect [g,f ] =1} : 
(w [f , t , sec_ord [f , 1] ] - 
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w[g,t-turn[g] ,sec_ord[g,maxi [g]]] )<=0; 
^turnaround constraint 
/*defn of arrived by*/ 

s.t. arrvby {f in F,k in K,t in T: f eas [t ,k,f ] =1 }: 

(w[f ,t,k] -w[f ,t-l,k])>=0; 

end; 

6.3 Sample data for the model 

The data file can be maintained separately from the model, to permit different data 
to be solved using the same model. The data has to be in a certain format for the 
interpreter-solver (glpsol) to work properly. Refer to [6] for details. 

6.3.1 The data file (unformatted) 
data; 

# Mumbai Bangalore Pune Goa Calicut Belgaum Co 
imbatore ; 

# MPG 1 GC 2 PbB 3 GCB 4 

set K := Mu Ba Pu Go Ca Be Co; 

set F :=Mu_Pu_Go Go_Ca Pu_Be_Ba Go_Co_Ba; 

set T := 1 2 3 4 5 6 7 8 9 10; 

set N:= 1 2 3; 

param maxi : = 
Mu_Pu_Go 3 
Go_Ca 2 
Pu_Be_Ba 3 
Go_Co_Ba 3; 

param connect default 0.0 := 
[Mu_Pu_Go,Go_Ca] 1; 

param feas default 0.0 := 

[*,*,Mu_Pu_Go] :Mu Ba Pu Go Ca Be Co := 

1.11... 

1 1 .11... 

2 1. 11... 

3 1.11. . . 

4 . . 1 1 . . . 

5 . . . 1 . . . 

6 

7 
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8 
9 

10 



[*,*,Pu_Be_Ba] :Mu Ba Pu Go Ca Be Co := 

. . . 

1 . . . 

2 . . 1 

3 . . 1 . . 1 

4 . . 1 . . 1 

5 . 1 1 . . 1 

6 . 1 ... 1 

7 . 1 ... 1 

8 . . . 

9 . . . 

10 . . 



[*,*,Go_Ca] : Mu Ba Pu Go Ca Be Co := 

. 

1 . 

2 . 

3 ... 1 . 

4 ... 1 1 

5 ... 1 1 

6 ... 1 1 

7 .... 1 

8 .... 1 

9 . 
10 



[*,*,Go_Co_Ba] :Mu Ba Pu Go Ca Be Co := 



1 

2 . . 

3 . . . 1 . . . 

4 ... 1 .. 1 

5 . 1 . 1 . . 1 

6 . 1 .... 1 

7.1 

8.1 

9 

10 ; 
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param sec_ord := 

[Mu_Pu_Go,l] Mu 

[Mu_Pu_Go,2] Pu 

[Mu_Pu_Go,3] Go 
[Go_Ca,l] Go 
[Go_Ca,2] Ca 

[Pu_Be_Ba,l] Pu 

[Pu_Be_Ba,2] Be 

[Pu_Be_Ba,3] Ba 

[Go_Co_Ba,l] Go 

[Go_Co_Ba,2] Co 

[Go_Co_Ba,3] Ba; 



param arr default 2.0 := 



param dep default 2.0 := 



param cap default 2.0 := 



param tr: Mu Ba Pu Go Ca Be Co 



Mu_Pu_Go 1111111 
Go_Ca 1111111 
Pu_Be_Ba 1111111 
Go_Co_Ba 111111 1; 

param turn [Mu_Pu_Go] :=1; 

param sched_dep:= 

Mu_Pu_Go 1 
Go_Ca 4 
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Pu_Be_Ba 3 
Go_Co_Ba 3; 

param sched_arr:= 
Mu_Pu_Go 3 
Go_Ca 5 
Pu_Be_Ba 5 
Go_Co_Ba 5; 

param air := 
Mu_Pu_Go 20000 
Go_Ca 30000 
Pu_Be_Ba 40000 
Go_Co_Ba 50000; 

param ground := 
Mu_Pu_Go 400 
Go_Ca 600 
Pu_Be_Ba 800 
Go_Co_Ba 1000; 

param atleast_one default 0.0 := 

[Mu_Pu_Go , * , *] Mu 3 1 

Pu 4 1 

Go 5 1 
[Go_Ca,*,*] Go 6 1 
Ca 7 1 

[Pu_Be_Ba,*,*] Pu 4 1 
Be 6 1 
Ba 7 1 

[Go_Co_Ba,*,*] Go 4 1 
Co 6 1 
Ba 7 1; 

end; 

7 Solution 

The LP relaxation of the TFMP model is first solved with Vwy 4 G 3? PI [0, 1] to verify 
if our data set supports the claims made in the paper by Bertsimas et.al [2]. 
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7.1 Solution to the LP relaxation of TFMP with conflicting 
flights 

A set of imaginary flights Pu_Be_Ba and Go_Co_Ba having the same scheduled arrival 
timings at Bangalore, which has an arrival capacity for only one flight was created, 
as in Table 1. However, even for this conflicting data set, the LP relaxation still 
had integral solutions, lending support to the results due to Bertsimas et.al. The 
Pu_Be_Ba flight is held for one interval after its scheduled departure from Goa, in 
order to satisfy the conflicting constraints. The minimal cost is 800 corresponding to 
the value of the Pu_Be_Ba flight, c 9 Pu Be Ba . 



Flight 


Departure time 


Arrival time 


Ground delay 


Air delay 


Scheduled 


Actual 


Scheduled 


Actual 


Mu_Pu_Go 


1 


1 


3 


3 








Go.Ca 


4 


4 


5 


5 








Pu_Be_Ba 


3 


4 


5* 


6 


1 





Go_Co_Ba 


3 


3 


5* 


5 









All airports have arrival capacity =1, at all times 



Table 1: Solution to the LP relaxation of TFMP with conflicting flights. 

Conflicting data set with PuJBeJBa and Go.Co.Ba having the same scheduled arrival timings 
(as shown by the asteriks (*)) at Bangalore, which can accomodate only one arrival at all 
time intervals. Notice in the following solution, how the Pu-BeJBa flight suffers a ground 
delay of 1 unit to resolve the conflict. 



Problem: lp_tfmp 

Rows : 588 

Columns: 316 

Non-zeros : 582 

Status : OPTIMAL 

Objective: cost = 800 (MINimum) 

No. Row name St Activity Lower bound Upper bound Marginal 



1 cost 



800 



No. Column name St Activity Lower bound Upper bound Marginal 



1 w[Mu_Pu_Go,0,Mu] 
B 
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2 w[Mu_Pu_Go,l,Mu] 

B 

3 w[Mu_Pu_Go,2,Mu] 

B 

4 w[Mu_Pu_Go,3,Mu] 

B 

5 w[Go_Ca,2,Go] 

B 

6 w[Go_Ca,3,Go] 

B 

7 w[Go_Ca,4,Go] 

B 

8 w[Go_Ca,5,Go] 

309 grdel_fl[Mu_Pu_Go] 

B 

310 grdel_fl[Go_Ca] 

B 

311 grdel_fl[Pu_Be_Ba] 

B 

312 grdel_fl[Go_Co_Ba] 

B 

313 airdel.f 1 [Mu_Pu_Go] 

NL 

314 airdel.f 1 [Go.Ca] 

B 

315 airdel.f 1 [Pu_Be_Ba] 

B 

316 airdel.f 1 [Go_Co_Ba] 

B 



1 
1 
1 



1 



1 


























19600 



Karush-Kuhn-Tucker optimality conditions: 



KKT . PE : max . abs . err . 

max . rel . err . 
High quality 



0.00e+00 on row 
0.00e+00 on row 



KKT . PB : max . abs . err . 

max . rel . err . 
High quality 



0.00e+00 on row 
0.00e+00 on row 



KKT.DE: max. abs. err. = 0.00e+00 on column 
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max.rel.err. = 0.00e+00 on column 
High quality 

KKT.DB: max.abs.err. = 0.00e+00 on row 
max.rel.err. = 0.00e+00 on row 
High quality 

End of output 

7.2 Baseline solution without conflicts 

The LP relaxation of the TFMP is now solved with the baseline data set, as in Table 2. 
As can be seen from the solution below, the linear program adjusts the departure and 
arrival timings so that the total cost is minimum. For the baseline data set, the cost 
is zero, implying that there are no delays of any sort. 



Flight 


Departure time 


Arrival time 


Ground delay 


Air delay 


Scheduled 


Actual 


Scheduled 


Actual 


Mu_Pu_Go 


1 


1 


3 


3 








Go_Ca 


4 


4 


5 


5 








Pu_Be_Ba 


3 


3 


5* 


5 








Go_Co_Ba 


3 


3 


5* 


5 









All airports have arrival capacity =2, at all times 



Table 2: Baseline solution without conflicts. 

Note the baseline data set is conflict free, since though Pu_BeJBa and Go-CoJBa have the 
same scheduled arrival timings (as shown by the asteriks (*)) at Bangalore, both can be 
accomodated due to the increased arrival capacity of 2. Notice that all flights are delay-free, 
in this instance. 



Problem: tfmp 

Rows : 588 

Columns: 316 

Non-zeros : 582 

Status : OPTIMAL 

Objective: cost = (MINimum) 

No. Row name St Activity Lower bound Upper bound Marginal 

1 cost B 
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No. Column name St Activity Lower bound Upper bound Marginal 



1 w[Mu_Pu_Go,0,Mu] 

B 

2 w[Mu_Pu_Go, l,Mu] 

B 

3 w[Mu_Pu_Go,2,Mu] 

B 

4 w[Mu_Pu_Go,3,Mu] 

B 

5 w[Go_Ca,2,Go] 

B 

6 w[Go_Ca,3,Go] 

B 

7 w[Go_Ca,4,Go] 

B 

8 w[Go_Ca,5,Go] 

NU 

309 grdel_fl[Mu_Pu_Go] 

B 

310 grdel_fl[Go_Ca] 

B 

311 grdel_fl[Pu_Be_Ba] 

B 

312 grdel_fl[Go_Co_Ba] 

B 

313 airdel.f 1 [Mu_Pu_Go] 

NL 

314 airdel.f 1 [Go.Ca] 

B 

315 airdel.f 1 [Pu_Be_Ba] 

B 

316 airdel.f 1 [Go_Co_Ba] 

B 




1 
1 
1 


1 
1 
































-600 



19600 



Karush-Kuhn-Tucker optimality conditions: 

KKT.PE: max.abs.err. = 0.00e+00 on row 
max.rel.err. = 0.00e+00 on row 
High quality 
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KKT.PB: max . abs . err . = 0.00e+00 on row 
max . rel . err . = 0.00e+00 on row 
High quality 

KKT.DE: max. abs. err. = 0.00e+00 on column 
max. rel. err. = 0.00e+00 on column 
High quality 

KKT.DB: max. abs. err. = 0.00e+00 on row 
max. rel. err. = 0.00e+00 on row 
High quality 

End of output 

7.3 Arrivals/Departures before scheduled time: negative costs 

A further interesting instance is when the actual arrival/departure timings are before 
the scheduled timings. Clearly, there is nothing restrictive about the formulation 
which will prevent such an instance from being considered. However, the economic 
interpretation is not the same. Moreover, with the objective being of minimization, 
the solver forces high values of negative air delays. In fact, there might be a situation, 
where one flight has a very high negative air delay and others have low positive values 
of ground delay. The real total cost, then is not negative. Actually, wherever a cost 
variable has a negative value it can be substituted by an auxiliary variable indicating 
its real value as 

a > max(0, 

where Real cost= a and computed cost=/3. The objective function will be (3 while 
the problem is to minimize it. The solution sumarized in Table 3, indicates such an 
instance of a negative cost, due to a negative air delay for the Go_Ca__Ba flight. 

Problem: tfmp 

Rows: 591 

Columns: 318 

Non-zeros : 594 

Status : OPTIMAL 

Objective: dummy = -2000 (MINimum) 



No. Column name St Activity Lower bound Upper bound Marginal 



1 w[Mu_Pu_Go,0,Mu] 
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B 

2 w[Mu_Pu_Go,l,Mu] 

B 1 

3 w[Mu_Pu_Go,2,Mu] 

NU 1 

4 w[Mu_Pu_Go,3,Mu] 

NU 1 

5 w[Go_Ca,3,Go] 

B 
309 grdel_fl[Mu_Pu_Go] 

B 

310 grdel_fl[Go_Ca] 

B 

311 grdel_fl[Pu_Be_Ba] 

B 

312 grdel_fl[Go_Co_Ba] 

B -2 

313 airdel.f 1 [Mu_Pu_Go] 

B 

314 airdel.f 1 [Go_Ca] 

B 

315 airdel.f 1 [Pu_Be_Ba] 

B 

316 airdel.f 1 [Go_Co_Ba] 

B 

317 cost B 











-400 
-58800 



Flight 


Departure time 


Arrival time 


Ground delay 


Air delay 


Scheduled 


Actual 


Scheduled Actual 


Mu_Pu_Go 


1 


1 


3 3 








Go_Ca 


4 


4 


5 5 








Pu_Be_Ba 


3 


3 


5 5 








Go_Co_Ba 


3 


1* 


5 3 


-2 






All airports have arrival capacity =2, at all times 



Table 3: Arrivals/Departures before scheduled time: negative costs. 
The feasible set of times Tq _Co_Ba is extended backward in time by 2 intervals, for all 
airports and sectors for the Go_Co_Ba flight. This permits the flight to depart before time 
from Goa airport, (see asterik (*)) and enter all sectors before its scheduled times, finally 
arriving before-time at Bangalore airport. In particular, note how the entire freedom-to- 
adjust (of 2 intervals) available in the schedule is taken up by the solution. 
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318 dummy obj 



-2000 



Karush-Kuhn-Tucker optimality conditions: 



KKT . PE : max . abs . err . 

max . rel . err . 
High quality 



5.68e-12 on row 4 
2.27e-13 on row 2 



KKT . PB : max . abs . err . 

max. rel . err . 
High quality 



4.44e-16 on column 14 
2.22e-16 on row 579 



KKT . DE : max . abs . err . 

max. rel . err . 
High quality 



6.26e-10 on column 16 
3.71e-10 on column 305 



KKT . DB : max . abs . err . 

max . rel . err . 
High quality 

End of output 



0.00e+00 on row 
0.00e+00 on row 



7.4 Departure of an outgoing flight scheduled before arrival 
of the incoming connecting flight 

This solution, seen in Table 4, indicates how the Go_Ca flight is delayed on the ground 
at Goa airport, because its scheduled departure is before the scheduled arrival of its 
incoming-connecting Mu_Pu_Go flight at Goa airport. 

Problem: tf 

Rows : 585 

Columns: 316 

Non-zeros: 560 

Status : OPTIMAL 

Objective: cost = 1200 (MINimum) 



No. Column name St Activity Lower bound Upper bound Marginal 

1 w[Mu_Pu_Go,0,Mu] 

B 1 

2 w[Mu_Pu_Go, l,Mu] 
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B 

3 w[Mu_Pu_Go,2,Mu] 

B 

4 w[Mu_Pu_Go,3,Mu] 

B 

5 w[Go_Ca,3,Go] 

B 

6 w[Go_Ca,4,Go] 

B 

7 w[Go_Ca,5,Go] 

B 

8 w[Go_Ca,6,Go] 

B 

9 w[Pu_Be_Ba,2,Pu] 

B 

10 w[Pu_Be_Ba,3,Pu] 
B 

309 grdel_fl[Mu_Pu_Go] 

NL 

310 grdel_fl[Go_Ca] 

B 

311 grdel_fl[Pu_Be_Ba] 



1 
1 
1 



0.5 
0.5 
1 

1 

2 

















1000 



Flight 


Departure time 


Arrival time 


Ground delay 


Air delay 


Scheduled 


Actual 


Scheduled Actual 


Mu_Pu_Go 


1 


1 


3* 3 








Go_Ca 


2* 


4 


4 6 


2 





Pu_Be_Ba 


3 


3 


5 5 








Go_Co_Ba 


3 


3 


5 5 









All airports have arrival capacity =2, at all times 



Table 4: Departure of an outgoing flight scheduled before arrival of the incoming 
connecting flight. 

Here, at Goa airport, the departing flight Go-Ca, is scheduled to depart, before the scheduled 
arrival of its incoming, connecting flight Mu.Pu.Go, arriving from Mumbai. The asteriks 
(*) indicate a situation where passengers getting off the MuJPu_Go flight at Goa airport 
would miss their connecting flight, Go_Ca, to Calicut. Notice how the LP resolves this 
issue by holding the Go-Ca flight on the ground for 2 time intervals, beyond its scheduled 
departure, so as to accomodate passengers from the incoming flight. Note that the schedule 
and connectivity parameters are such, that this delay is inevitable, and the LP ensures that 
exactly those flights are delayed which minimize the cost of delay. 
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B 

312 grdel_fl[Go_Co_Ba] 

B 

313 airdel.f 1 [Mu_Pu_Go] 

NL 

314 airdel.f l[Go_Ca] 

NL 

315 airdel.f 1 [Pu_Be_Ba] 

B 

316 airdel.f 1 [Go_Co_Ba] 

B 





















20600 
30000 



Karush-Kuhn-Tucker optimality conditions: 



KKT . PE : 



max . abs . err . 
max . rel . err . 
High quality 



0.00e+00 on row 
0.00e+00 on row 



KKT . PB : max . abs . err . 

max . rel . err . 
High quality 



0.00e+00 on row 
0.00e+00 on row 



KKT . DE : max . abs . err . 

max . rel . err . 
High quality 



0.00e+00 on column 
0.00e+00 on column 



KKT.DB: max. abs. err. = 0.00e+00 on row 
max. rel. err. = 0.00e+00 on row 
High quality 



End of output 



8 Conclusions 

We have tried to explore the logic behind the claims by Bertsimas et.al about integral 
solutions to the LP relaxation of the TFMP. 

Polyhedral theory only indicates that the stronger TFMP formulation of Bertsimas 
et.al might lead to integral solutions in some cases. Unless we obtain all facet-defining 
inequalities of the conv(S), we are not assured of obtaining integral solutions from 
the LP relaxation. Therefore, only computations can provide a rough estimate of the 
frequency of integral solutions, in the worst case, which in turn will decide whether 
the formulation offers practical advantages. Our computations indicate that the en- 
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couraging results reported by Bertsimas et.al are not merely fortuitous or due to their 
specific data set. Indeed, we found that the TFMP had integral solutions even in case 
of artificial data sets generated to include severe conflicts in the flight schedules. In 
our limited tests with 4-5 scenarios, we obtained non-integral solutions only once. 

This is of significant practical importance because, the LP relaxation can be solved 
even on small machines with low memory and processor speed. The formulation with 
308 variables took less than 0.1 seconds when solved on a regular IBM laptop having 
1GB RAM and 1.66 GHz Processor speed. In contrast, the Integer program has to be 
solved by Branch and Bound which would be much more expensive and difficult to 
implement on such a machine, especially since the number of variables rises rapidly 
with addition of new flights. The importance of obtaining integral solutions, using the 
LP relaxation itself (without the need for Branch and Bound), will be more acutely 
felt for problems of large sizes. 

9 Directions for future work 

An important observation, is that the size of the problem is quite large even for simple 
cases. E.g. for a data set with 4 flights, 7 sectors and 11 time slots, the formulation 
already has 308 variables. In a real life data set, the proportion of conflicting flights 
may not be very high. In such a situation, it is better to solve the LP for only those 
flights which are in conflict with each other. 

One approach could be, to use the scheduled arrival and departure timings,to 
generate a set of flights which are in conflict with each other, called the set SC . The set 
W = &\3£ then consists of flights which are not in conflict either among themselves 
or with those in S£ . The TFMP is solved only for flights / G SC . Then, the entire set 
& is again checked for conflicts, with the flight timings for those in 3£ now updated 
to those given by the TFMP solution. If this causes some flights in & to be in conflict 
with those in 2C , such flights are transferred from <¥ to 3£ . This iterative process 
is continued till we get a schedule with no conflicts. 

The problem with such an approach is that the set S£ might grow, until finally 
2£ = & thus offering no advantages in terms of problem size. Moreover, the TFMP 
has to be solved once per iteration and further analysis needs to be done to weigh 
the benefits of solving a reduced problem, against the expense in solving the TFMP 
several times. The concept, is however very attractive for an initial schedule having 
very few conflicting flights, if some heuristics can be developed to locally adjust the 
conflicting set 5£ without affecting the larger, non-conflicting set W . 
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